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Abstract 

We show that the occurrence of approximate pseudo-spin symmetry in nu- 
clei is connected with certain similarities in the relativistic single-nucleon wave 
functions of the corresponding pseudo-spin doublets. We perform a case study 
in which several examples and the systematics of this connection are explored. 
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I. INTRODUCTION 



In a recent paper [|l|] it was shown that quasi-degenerate pseudo-spin doublets in nu- 
clei, discovered almost thirty years ago arise from the near equality in magnitude 
of attractive scalar, Vs, and repulsive vector, Vy, relativistic mean fields, Vs ~ — Vy, in 
which the nucleons move. Pseudo-spin doublets have non- relativistic quantum numbers (n r , 
£j — ^ + 1/2) and (n r — 1, £ + 2, j = £ + 3/2) where n r ,£, and J are the single-nucleon 
radial, orbital, and total angular momentum quantum numbers, respectively. These authors 
defined a "pseudo" orbital angular momentum £ = £ + 1; for example, (n r sy2, (n r — 1)^3/2) 
will have £ — 1 , (n r p 3 / 2 , (n r — 1)^5/2) will have I — 2, etc. Then these doublets are almost 
degenerate with respect to "pseudo" spin, s = 1/2, since j = £ ± s for the two states in 
the doublet. This symmetry has been used to explain a number of phenomena in nuclear 
structure 0] including most recently the identical rotational bands observed in nuclei ||, 
and to establish an effective shell model coupling scheme |J. A near equality in the mag- 
nitude of mean fields seems to be a universal feature of relativistic theories ranging from 
relativistic field theories with interacting nucleons and mesons 0, to nucleons interacting 
via Skyrme-type interactions to QCD sum rules JTIJ. In this paper we investigate the 
consequences of Vs ~ —Vy on the single-nucleon wave functions. The connection between 
the Dirac equation and pseudo-spin symmetry is reviewed in Sec. II. The relativistic mean 
field model described in Sec. Ill is used for the calculations of the case study presented in 
Sec. IV. Our conclusions are given in Sec. V. 



II. THE DIRAC EQUATION AND PSEUDO-SPIN SYMMETRY 

The Dirac equation with external scalar, Vs, and vector, Vy, potentials is given by: 

[ca-p + f3(mc 2 + V s ) + V v }y = EV, (1) 



where a and (3 are the usual Dirac matrices |TTJ. We shall be considering relativistic 



mean field theory with spherical symmetry for which the scalar and vector potentials de- 



pend only on the radial coordinate. In this case the orbital angular momentum is not 
a conserved quantum number in general. Instead a nucleon moving in a spherical rel- 
ativistic field is labeled by a radial quantum number, n r , total angular momentum j, 
its projection on the z-axis, m, and k = —(3 (a ■ L + 1) [11]]. The eigenvalues of k are 
k = ±(j + 1/2); - for aligned spin (si/2,P3/2, etc.) and + for unaligned spin G»i/2,<^3/2, 
etc.). Thus, the quantum number k and the radial quantum number n r are sufficient to 
label the orbitals. The spherically symmetric Dirac wave function can then be written 
in terms of upper and lower components *«>o,m = (5 , K[^x]m =£ ~ 1/2) ^/K[^-iX]m =£ ~ 1/2) ) 5 
^K<o,m = (#4^x]m =m/2 \z/4^mX]m =m/2) ) where g K , f K are the radial wave functions 
(omitting the radial quantum numbers), Yi are the spherical harmonics, x is a two- 
component Pauli spinor, and [. . means coupled to angular momentum j. The radial 
Dirac equation for the upper and lower components (g K , f K ) of the single-nucleon radial 
wave function in dimensionless units (also omitting the magnetic quantum numbers) is then 
given by |TT 



d 1 + k 
dr r 



g K = [2 - E - V{r) }f K 



(2) 



d 1 — AC 
dr r 



f K = [E + A(r) ]g K , 



(3) 



where r is the radial coordinate in units of hc/mc 2 , V(r) = [Vy(r) - l^(r)] / mc 2 , A(r) 
= [Vg(r) + VV(r)] / mc 2 , and E is the binding energy [E = (mc 2 — S)/mc 2 > 0] of the 
nucleon in units of the free nucleon mass. First we repeat the proof [|IJ that, in the limit 
of equality of the magnitude of the vector and scalar potential, A(r) = 0, pseudo-spin is 
exactly conserved. To do this, we solve for g K in (|3|) and substitute into (H), obtaining the 
second order differential equation for f K , 



d 2 



2 d 



dx 2 x dx 



+ 1) 



x. 



+ (V(r) — 2 + E) 



/« = 0, 



(4) 



where x = \[E r, E 7^ 0, and 



K — 1, K > ; i — —K, K < 0, 



(5) 



which agrees with the original definition of the pseudo-orbital angular momentum 0,[3] . For 
example, for (n r Si/ 2 , (n r — 1)^3/2), K — — 1 and 2, respectively giving £ = 1 in both cases. 
Furthermore, the physical significance of £ is revealed; it is the "orbital angular momentum" 
of the lower component of the Dirac wave function and, in this limit, it is a conserved 
quantum number. 

For E 7^ 0, Equation (Q) is a Schrodinger equation with an attractive potential V and 
binding energy 2 — E which depends only on the pseudo-orbital angular momentum, £, 



through the pseudo - rotational kinetic energy, 



„2 , 



and not on n. Hence the eigenenergies 



do not depend on k but only on £. Thus the doublets with the same £ but different k (k, = £+1 
and n = —£) will be degenerate, producing pseudo-spin symmetry. 

Furthermore, from it is clear that the f K with the same £ can only differ by an 
overall constant, 



f-i = Ml f\ 



(6) 



We show now that Mf = 1. The Dirac wave function is normalized by 



r 2 dr[g K (r) 2 + f K (rf] = 1. 



Furthermore, from (H), we find in the pseudo- spin limit, 



(7) 



1 

Ei 



d £ 
dr r 



1 



d £ + 1 
dr r 



f-l(r) 



(8) 



where we have used E~ t = Ei +l = E_~ t . Using these relations in (^) for both partners, k = — £ 
and k = £ + 1 and integrating by parts and then using the differential equation (^) for f K , 
we obtain 

l r 2 drf l+1 (rY(2-V-2E l ) = -E i , 

Mf / r 2 drf l+1 (r) 2 (2 - V - 2E~) = -E b (9) 

which implies that j\T~ = 1 for E^ 7^ 0. Furthermore, for a square well |TJ, A/| = — 1 and 
therefore we assume that this is the case for more realistic nuclear mean fields as well. 



As pointed out in [lj], there will not be any bound Dirac valence states, only Dirac 
sea states, for the pseudo-spin limit A = 0, which would contradict the fact that bound 
nuclei exist. In |I| it was shown that it is possible to have bound valence states and quasi- 
degeneracy for a small A(r) by studying the Coulomb potential and the spherical potential 
well, = V s ,v > , r < R ; V s>v (r) = , r > R. Likewise, one may ask how well the 

relation 

Li(r) = -f l+1 (r) (10) 

holds for small pseudo-spin breaking. In this paper we employ a physically more realistic 
Dirac-Hartree approach to answer this question. 



III. RELATIVISTIC POINT COUPLING MODEL 

We use a self-consistent Dirac-Hartree model with contact interactions (point couplings) 
in the mean field (ip — > (V ; )) an d no Dirac sea approximations. The model consists 
of four-, six-, and eight-fermion point couplings leading to scalar and vector densities with 
both isoscalar and isovector components, derivatives of the densities to simulate the finite 
ranges of the mesonic interactions, but no explicit meson fields. The Lagrangian is given by 

C — Cf ree + £if + Chot + £der + £em , where (11) 
Cjree and C em are the kinetic and electromagnetic terms, and 

£ 4/ = -~as(W>)(iH>) ~ g ^(^7^0(^7^) 
1 _ _ l 



Chat = -\f3s(^f ~ \ls(^Y 

-\lvBl^m 2 , and (13) 
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-\s v {d v i>i^W4>ri>) • (14) 

In these equations, ip is the nucleon field, the subscripts "S"' and "V" refer to the scalar 
and vector nucleon fields, respectively, and the subscript "T" refers to isovector fields con- 
taining the nucleon isospin f. The physical makeup of C is that £4/ is a four-fermion inter- 
action, while Chot contains higher order six-fermion and eight-fermion interactions, and Cd er 
contains derivatives in the nucleon densities. There are a total of nine coupling constants. 

Minimizing the expectation value of the Hamiltonian corresponding to Eq. (11) in the 
space of Slater determinants \(f>) leads to the Dirac-Hartree equations containing the following 
potentials: 



V s = asps + PsPs + lsP 3 s + Ss^ps , (15) 

V v = a v p v + -fvpy + 5 v Ap v , (16) 

V T s = oitsPts , and (17) 

Vtv = OL TV p T y , (18) 



where Eq. (15) is the isoscalar-scalar potential corresponding to a meson (fictitious) ex- 
change, Eq. (16) is the isoscalar-vector potential corresponding to uj meson exchange, Eq. 
(17) is the isovector-scalar potential corresponding to 5 meson exchange, and Eq. (18) is the 
isovector-vector potential corresponding to p meson exchange. In these latter equations the 
scalar density is given by ps = the vector density is given by py = (0|^7o^|</>), 

the isovector-scalar density is given by pxs = (0|'07'3'0|0)> and the isovector-vector density 
is given by prv = (01 ^7370^ !</>)• These densities can of course be rewritten in terms of 
the upper and lower components of the Dirac wave functions (g a ,fa), respectively, where 
a = {n r , K,m}. For example, the vector (baryon) density becomes 

occ - Q-j'-l-l 

Mr) = E J £-[^(0 + ^(01, (19) 

while the scalar density becomes 
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° cc - On _|_ 1 

Ps(r)=J2^-[gl(r)-f a (r)}- (20) 

The nine coupling constants of the Lagrangian were determined in a self-consistent pro- 
cedure that solves the Dirac-Hartree equations for several nuclei simultaneously in a nonlin- 
ear least-squares adjustment algorithm of Levenberg-Marquardt type with respect to well- 
measured nuclear ground-state observables. The well-measured observables used (in order 
of decreasing experimental accuracy) are (a) the ground-state masses (binding energies), (b) 
the rms charge radii, and (c) the spin-orbit splittings of the least-bound neutron and proton 
spin-orbit pairs. We chose these observables for ten spherical closed-shell nuclei ranging from 
16 to 208 Pb in the determination of the coupling constants (40 observables to determine 9 
coupling constants). The final coupling constants |9j are given in Table I where the first four 
refer to Eq. (12), the next three refer to Eq. (13), and the remaining two refer to Eq. (14). 
They span 13 orders of magnitude, but seven of them and the sum of the remaining two are 
natural (dimensionless numbers of order 1) when scaled in accordance with QCD mass scales 
and taking into account the constraint of chiral symmetry |9|,0. They should be viewed 



as an interim set of coupling constants in the development of the relativistic point coupling 
model. Note that the initial set of coupling constants determined with this approach appears 
in Ref. §. 

With these nine coupling constants one can calculate the following for spherical closed- 
shell nuclei: (a) Dirac single-nucleon wave functions and eigenvalues for both protons and 
neutrons, (b) nuclear ground-state mass and binding energy, (c) proton and neutron densities 
and their moments, (d) nuclear charge density and its moments, and (e) isoscalar- and 
isovector-, scalar and vector, potentials. 

IV. CASE STUDY 

For our present purpose we consider the upper and lower components of the Dirac single- 
nucleon wave functions for 208 Pb. Using the coupling constants of Table I we have calculated 
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the single-neutron wave functions for the occupied orbits in the ground state of this nucleus; 
these orbits have six neutron pseudo-spin doublets. Their properties are listed in Table II. 
Similarly, we have calculated the single-proton wave functions for the occupied orbits in the 
ground state of the same nucleus which contain four proton pseudo-spin doublets, and their 
properties are listed in Table III. The pseudo-spin doublets are ordered in the two tables 
according to increasing orbital angular momentum £ of the lower Dirac component f nrjK , 
for each value of the radial quantum number n r . The energy eigenvalues therein, e nrtK , in 
units of MeV, are related to the dimensionless energy eigenvalues E HrjK , by 

£n r ,n = E nrtK mc 2 (21) 

where mc 2 is either the free neutron mass or the free proton mass. 

The upper and lower components of the Dirac wave functions for the most deeply bound 
neutron and proton pseudo-spin doublets (first entries in Tables II and III, respectively) are 
shown in Figs. 1-4. In Figs. 2 and 4 the lower components of the doublets are plotted 
with opposite sign so as to directly test Eq. (10), the pseudo-spin limit, against the observed 
(calculated) pseudo-spin breaking. It is clear that Eq. (10) is crudely satisfied when the 
eigenvalues are crudely degenerate (to within ~ 3 MeV for both neutrons and protons in 
this example). Comparisons of Figs. 1 and 2, for neutrons, and Figs. 3 and 4, for protons, 
show that while the upper components of the total Dirac wave function dominate, the lower 
components are certainly not negligible, and approximate, although crude in this example, 
pseudo-spin symmetry exists for both protons and neutrons. 

We examine the systematics of pseudo-spin symmetry and the corresponding relativistic 
wave functions by considering the six neutron pseudo-spin doublets in Table II. Figures 5 - 
7 show the lower components of three of the remaining five doublets. In Figs. 5 and 6, for 
£ = 2 and 4, respectively, the pseudo-spin limit given by Eq. (10) is again approximately 
satisfied to about the same degree as that for the first neutron pseudo-spin doublet for £ 
= 1 shown in Fig. 2. The same remark applies to the case for £ = 3 (not shown). The 
four doublets all have the same radial quantum number n r — 1, and have increasing orbital 
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angular momentum £ of the lower component. Figure 7 shows the lower components for the 
first of the two neutron pseudo-spin doublets that have radial quantum number n r = 2. Here 
there is dramatic improvement in pseudo-spin symmetry, as expected [|| because these states 
are less bound and exact pseudo-spin symmetry occurs for no bound Dirac valence states. 
The same is true for the second neutron pseudo-spin doublet with n r = 2 in Table II (not 
shown). A related systematic is that if the Dirac single- nucleon states of the pseudo-spin 
doublets of Tables II and III are sorted by classical major oscillator shells, then for those 
states in the same major shell, the energy splitting of the doublets decreases as £ decreases. 
In Table II, the (n r ,£) = (1,3) and (2,1) doublets are in the N = 4 major shell and the 
(1,4) and (2,2) doublets are in the N = 5 major shell. In both cases the energy splitting 
decreases with decreasing i. The same is true for the (1,3) and (2,1) doublets in Table III 
also corresponding to the N = 4 major shell. This point becomes more clear in Fig. 8 which 
shows the spectra for the twelve states making up the six neutron doublets and the eight 
states making up the four proton doublets. Clearly, for the doublets with the same pseudo- 
orbital angular momentum, the splitting between doublets decreases as the binding energy 
decreases, or, equivalently, as the radial quantum number increases. Also for doublets with 
different pseudo-orbital angular momentum, but roughly the same binding energy, say for 
example, (n r ,£) = (1,3) and (2,1), and (n r ,£) = (1,4) and (2,2), the doublet with the lower 
£ has the smaller energy splitting. These features are seen in the square well potential as in 
Equation (19) in in which the energy splitting between the doublets is proportional to 
2£ + 1 and E. The six neutron pseudo-spin doublet splittings of Fig. 8 are again plotted 
in Fig. 9, with the factor (2£ + 1) divided out, vs the mean eigenvalue < e > for each 
doublet. Here it is clearly quantitative that the energy splitting decreases as the binding 
energy decreases. 

Finally, we compare two calculated neutron pseudo-spin doublets with experiment in 
Table II and two calculated proton pseudo-spin doublets with experiment in Table III. For 
both protons and neutrons, the magnitudes of the calculated eigenvalues and corresponding 
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measured values are comparable and the calculated energy splittings agree to within a factor 
~2 with the measured energy splittings in each of the four cases. The calculated splittings 
are always larger than the measured splittings which is not very surprising as most bound- 
state single-nucleon eigenvalue calculations, in relativistic mean field approaches, share the 
common feature of "spread-out eigenvalues" (see, for example, Ref. and our interim 

coupling constants (Table I) provide no exception. In addition, the measured and calculated 
states making up the neutron (n r ,£) = (2,2) pseudo-spin doublet in Table II are reversed, a 
feature also occurring in other relativistic calculations [[7U§ . On the other hand, the measured 
and calculated states making up the neutron (1,4) pseudo-spin doublet in the same table are 
not reversed in our calculation, but are reversed in other recent non-relativistic calculations 



T3|]l4}| . We remark here that these latter features do not change our observations as to 



the connections between (approximate) pseudo-spin symmetry and the corresponding Dirac 
single-nucleon wave functions and eigenvalues. 



V. CONCLUSIONS 

We have shown that physically realistic relativistic mean fields lead to small pseudo-spin 
symmetry breaking which implies that the lower components of the corresponding Dirac 
single-nucleon wave functions satisfy 

/_|(r) c- -f i+1 (r) (22) 

where the pseudo-orbital angular momentum £ is exactly conserved in the pseudo-spin limit. 
Near this limit the lower components depend only upon I (for fixed radial quantum number 
n r ) and the energy splitting between the doublets is small. Furthermore the pseudo-spin 
symmetry becomes increasingly valid as the pseudo-orbital angular momentum I decreases 
and as the binding energy Eg decreases (or as the radial quantum number n r increases). 
These observations and conclusions are based upon our calculations of ten pseudo-spin dou- 
blets occurring in 208 Pb of which four have been measured and compare reasonably well 
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with the corresponding calculated doublets. Thus this work confirms the contention of Ref. 
|Ij] that pseudo-spin symmetry in nuclei arises largely from nucleons moving in relativistic 
mean fields with attractive isoscalar-scalar and repulsive isoscalar-vector components that 
are nearly equal in magnitude. 
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TABLES 

TABLE I. Optimized Coupling Constants for the Interim Relativistic Point Coupling Model 



Coupl 



ing Constant 


Magnitude 


Dimension 


a s 


-4.517xl0~ 4 


MeV" 2 




-2.168xl(T 5 


MeV" 2 


ay 


3.435 xlO~ 4 


MeV" 2 


OLTV 


5.365xl0" 5 


MeV" 2 


Ps 


1.137X10 -11 


MeV" 5 


IS 


5.731xl0- 17 


MeV" 8 


IV 


-4.423xl0~ 17 


MeV~ 8 


6s 


-4.282xl0~ 10 


MeV" 4 


S v 


-1.155xl0- 10 


MeV" 4 



TABLE II. Calculated and Measured Neutron Pseudo-Spin Doublets in 208 Pb 



I 


n r , 


K< 




^n r ,K<0 


n r — 1, k> 




£n r — 1,«>0 










(MeV) 






(MeV) 


1 


1 


-1 


(«l/2) 


42.139 


0, 2 


{d-i/2) 


45.331 


2 


1 


-2 


(PS/2) 


31.454 


0, 3 


if 5/2) 


35.490 


3 


1 


-3 


(4/2) 


20.981 


0, 4 


(97/2) 


24.741 


4 


1 


-4 


(/7/2) 


10.944 


0, 5 


(hg/2) 


13.519 


4 


1 


-4 


(/7/2) 


9.708 a 


0, 5 


(hg/2) 


10.781 a 


1 


2 


-1 




18.129 


1, 2 


(^3/2) 


19.002 


2 


2 


-2 


(P3/2) 


7.656 


1, 3 


(/5/2) 


8.353 


2 


2 


-2 


(P3/2) 


8.266 a 


1, 3 


(fs/z) 


7.938 a 



^Experimental value from Ref. [15]. 
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TABLE III. Calculated and Measured Proton Pseudo-Spin Doublets in Pb 



I 


n r , k< 




^n r ,K<0 


n r — 1, k> 


0? + 2,i + l) 


£n r — 1,«>0 








(MeV) 






(MeV) 


1 


1,-1 


(*l/2) 


32.047 


0, 2 


(4/2) 


35.830 


2 


1,-2 


(P3/2) 


21.800 


0, 3 


(/s/2) 


26.358 


3 


1,-3 


(4/2) 


11.597 


0,4 


(57/2) 


15.930 


3 


1,-3 


(4/2) 


9.696 a 


0, 4 


(57/2) 


11.487 a 


1 


2, -1 


(«i/a) 


8.416 


1, 2 


(4/2) 


9.663 


1 


2, -1 


(si/2) 


8.013 a 


1, 2 


(4/2) 


8.364 a 



a Experimental value from Ref. fl5|| . 
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FIGURES 

FIG. 1. Dirac upper component wave functions for the neutron pseudo-spin doublet (n r ,£) = 
(1,1) in 208 Pb and specified more completely in the first line of Table II. 

FIG. 2. Dirac lower component wave functions corresponding to Fig. 1. The sign of the 
component with n > has been reversed to test Eq. (10). 

FIG. 3. Dirac upper component wave functions for the proton pseudo-spin doublet (n r ,£) = 
(1,1) in 208 Pb and specified more completely in the first line of Table III. 

FIG. 4. Dirac lower component wave functions corresponding to Fig. 3. The sign of the 
component with k > has been reversed to test Eq. (10). 

FIG. 5. Dirac lower component wave functions for the neutron pseudo-spin doublet (n r ,£) 
= (1,2) in 208 Pb and specified more completely in the second line of Table II. The sign of the 
component with n > has been reversed to test Eq. (10). 

FIG. 6. The same as Fig. 5 except for the pseudo-spin doublet (n r ,£) = (1,4) and specified 
more completely in the fourth line of Table II. 

FIG. 7. The same as Fig. 5 except for the pseudo-spin doublet (n r ,£) = (2,1) and specified 
more completely in the sixth line of Table II. 

FIG. 8. Spectra of the six neutron and four proton pseudo-spin doublets occurring in 208 Pb. 
Note that the pseudo-spin symmetry is more valid for the higher value of the radial quantum 
number. 

FIG. 9. Dirac energy splittings of the six neutron pseudo-spin doublets shown in Fig. 8, 
normalized by (2£ + 1), vs the mean energy eigenvalue < e > of the doublet. The line is a linear 
fit to the four cases with n r = 1. 
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